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Abstract 

This chapter seeks to outline a few basic problems in quantum statistical physics 
where recent experimental advances from the atomic physics community offer the hope of 
dramatic progress. The focus is on nonequilibrium situations where the powerful concepts 
and methods of equilibrium statistical physics and "linear response" theory (for small 
deviations from equilibrium) are not applicable. The problems discussed here are chosen 
in part because they have a high degree of "universality" or generality across different 
microscopic situations, as the major challenge in nonequilibrium statistical physics, both 
quantum and classical, has been to find principles as general as the basic principles of 
equilibrium statistical physics or linear response. 



1. Introduction 

The initial engagement of the condensed matter community with ultracold atomic 
physics during the first half of the past decade was due to the apparent similarity between 
the newly-created atomic systems and familiar (electronic) counterparts in the solid state. 
The two most notable examples of this trend were the observation of the superfluid to 
Mott insulator transition in an optical lattice pT] , and the creation of superfluid states of 
paired fermions, the neutral analogs of s-wave superconductors |2]. 

One feature of the latter half of the decade was a growing appreciation for the differ- 
ences between these two settings for quantum many body physics. This chapter concerns 
one of the most prominent, intriguing, and conceptually far-reaching of these: the out of 
equilibrium nature of many experiments in ultracold physics. Mechanisms of equilibra- 
tion in solid state systems are typically fast - on the picosecond timescal^- with the 
resulting equilibrium states essentially immortal. In contrast, the lower temperatures 
and energy scales of ultracold gases means that the corresponding processes can be ob- 
served at the millisecond level, and non equilibrium physics is a fact of life during the 
few second lifespan of most ultracold systems. 



^Glassy systems are of course a notable exception. 
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Furthermore, there is good reason to beheve that the mechanisms of equihbration 
in the ultracold domain are distinct from those that are genericaUy important for elec- 
trons in the sohd state. There are no phonons that can transfer energy into the lattice, 
(typically) no impurities to allow momentum to dissipate, and no spin-orbit interaction 
to mediate spin relaxation. We find ourselves therefore in an enviable (if unfamiliar) 
position: the simple Hamiltonians that we write for the degrees of freedom of interest 
are, to a good approximation, all there is. To take a simple example, an applied Zeeman 
field will typically polarize the spins of an electron system, while in an atomic system 
it leads only to Larmor precession (the quadratic Zeeman effect can elicit interesting 
effects, however). 

The study of equilibration in isolated systems (though this term should be used 
carefully, see below) described by strictly Hamiltonian dynamics is of course as old as 
statistical mechanics itself, dating back at least to Boltzmann. It is the possibility of 
studying non equilibrium quantum phenomena in a very simple setting that has led to 
a resurgence of interest in this problem of fundamental physics [3 H] . Additionally, the 
way in which a system comes to equilibrium is believed to be sensitive to a number 
of factors that are already the focus of many experimental investigations in ultracold 
systems, namely dimensionality, disorder, integrability of the underlying dynamics, and 
of course the initial conditions. 

In trying to draw together some of these threads for this chapter, we are faced with the 
familiar problem: equilibrium systems are all alike, while every non-equilibrium system 
is out of equilibrium in its own wa}j^ Nevertheless, most phenomena of current interest 
can be associated with one or more of the above four aspects. Let us take as an example 
an experiment that has acquired an iconic status as an illustration of out of equilibrium 
behavior in an atomic gas. In the 'quantum Newton's cradle' of Kinoshita et al. 
the evolution of the momentum distribution of arrays of harmonically confined Id Bose 
gases was studied over many periods of oscillation, after initially splitting the gas into 
two counterpropagating clouds. Even though each atom undergoes many thousands of 
collisions during this time, the momentum distributions do not relax to the equilibrium 
distribution (determined to be Gaussian for thermal clounds of the same rms momentum) . 
This is to be contrasted with thermalization in a three-dimensional gas, which occurs 
after a very small number of collisions. 

A plausible argument for this behavior is easily made if one assumes the applicability 
of the Boltzmann equation in its naive form (which is by no means clear), with a colli- 
sion integral describing two-particle collisions. In a gas of particles of equal mass, such 
collisions do not lead to a change in the distribution function, as particles either retain 
or exchange their momenta. Within the Boltzmann picture, a change in the distribution 
function requires at least three-body collisions to be accounted for, as was done in a 
number of recent works |S1 [7] . According to the discussion of Ref. [5] , once the suppres- 
sion of three-body scattering due to intcrparticle repulsion is accounted for, the resulting 
damping rates are consistent with the long times observed by Kinoshita et al.. 

The example of the quantum Newton's cradle vividly illustrates the role played by 
reduced dimensionality. Additionally, the experimental conditions closely approximate 
one of the best-known integrable many-body systems: the Id Bose gas with (5-function 



■^With apologies to Tolstoy. 
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interactiorj^ In this integrable system, three-body conision are absent, or more correctly 
factorize into successive two-body cohisions. The three-body coUisions that will eventu- 
ally cause the system to equilibrate arise because the Id model is only an approximation 
to experimental reality. The potential confining atoms to a tube naturally has higher 
modes of transverse excitation, and it is virtual processes involving these modes that 
gives rise to the three-body interactions [6] . In Section |3] we will see another example of 
thermalization due to the breaking of integrability. 

Equilibration is also often deemed to be particularly precarious in disordered systems, 
which appear at first to be the polar opposites of integrable systems. The phenomenon 
of localization due to quantum mechanical interference has been studied constantly since 
Anderson's work [101 more than fifty years ago. Recently it was suggested [H, that 
an isolated disordered interacting system can undergo a finite temperature localization 
transition. In the low temperature (localized) phase the system is unable to come to in- 
ternal equilibrium. This fascinating idea would seem to have a natural home in ultracold 
systems and will be discussed in Section |4j 

We begin our survey with the last of the themes identified above: the choice of initial 
conditions. A natural way to initiate out of equilibrium dynamics is to abruptly change 
some system parameter. In the classical description of phase ordering one speaks of 
a quench if the variation of this parameter would cause the equilibrium system to pass 
through a phase transition. By analogy, abruptly crossing a quantum phase transition has 
become known as a quantum quench^ an alliterative coinage apparently due to Calabrese 
and Cardy [T2], though recent papers use the term to refer to any abrupt change in the 
system. Section [2] gives an introduction to some of the general theoretical issues with 
reference to a simple model, the 'A0^' field theory. 

As indicated above, this survey is necessarily incomplete and many interesting aspects 
of non-equilibrium behavior are missing. For example, we will have nothing to say about 
non-equilibrium steady states, which seem less relevant to ultracold physics than to 
(say) mesoscopic physics, where conditions of constant drive can be sustained for long 
periods. We will only have a few brief comments about connections between quantum 
non-equilibrium physics and quantum information concepts such as entanglement. It is 
hoped that readers will find the chapter a useful effort to motivate and explain some key 
current questions in non-equilibrium physics that are ripe for investigation. 

2. Quantum Quenches 

2. 1 . Introduction 

The study of phase ordering kinetics has an extremely long history in classical statisti- 
cal physics 1131 . The main theoretical problem is to understand how a new phase appears 
after an abrupt change in system parameters. For a symmetry-breaking transition, this 
entails an understanding of the growth of the order parameter, and the dynamics of topo- 
logical defects formed in the process. It is natural to ask what changes when we cross 
a quantum - rather than classical - phase transition, and more generally how quantum 
mechanics changes things. 



^Also known as the Lieb-Liniger gas |9]. 
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We would like to argue that for a symmetry-breaking transition, it is only in the 
early stages of phase ordering that the difference will be acutej^ The reason is that the 
initial dynamics establishes a local fluctuating order parameter distribution that soon 
becomes macroscopic in character and leaves the quantum domain. The subsequent 
evolution of this distribution (often called coarsening in the literature) will then be 
described by classical equations of motior[^ Thus the following discussion will focus on 
the quantum dynamics of the early stages of phase ordering and the emergence of the 
classical description within a simple model. We should point out that this issue was 
discussed many years ago in a cosmological context [13] , predating the more recent work 
motivated by the experiment of Ref. [15] ITTl ITSl [T9] . 

2.2. A model system 

Let try to describe the effect of a quench in a system described by the 'A^''' Hamil- 
tonian 



where $(r) and n(r) are complex canonically conjugate fields [$(r),n(r')] — iS{r — r'). 
This model describes a quantum phase transition between a disordered phase (($) = 0) 
at positive to an ordered one (($) ^ 0) at negative values. For simplicity, we imagine 
that initially the system is in the ground state corresponding to = mf > 0. 

Ultracold physics provides two examples of systems with transitions described by 
the model Eq. ([T]). The first is the superfluid insulator transition in the Bose-Hubbard 
model [20l [1], where a non-zero order parameter ($) corresponds to the condensate 
wavefunction in the superfiuid phase, and the transition may be tuned by varying the 
depth of the optical lattice that confines the bosons. The second example is provided by a 
condensate of spin-1 bosons with ferromagnetic spin-spin interactions. If the system has 
zero total magnetization with respect to the axis defined by the applied magnetic field, 
only the quadratic part of the Zeeman energy is effective and leads, at large fields, to a 
condensate in the m = state. As the field is lowered, the ground state of the condensate 
makes a transition to a state with nonzero transverse magnetization [51] [551 [IS] ■ If we 
consider only small changes of the magnetic field around the critical value, this problem 
maps onto the 'A0^' theory where ($) gives the (amplitude and phase of the) transverse 
magnetization |16) . The groundbreaking experiment described in Ref. |15j explored the 
dynamics following a quench in this system. 

Thus both realizations have a natural control parameter corresponding to m^, that 
can be quenched from positive to negative values. The 'A(/)^' theory is of course one of 
the first examples of a quantum field theory that one meets, and a prototype for studies 
of spontaneous symmetry breaking. Therefore it is not surprising that the dynamics of 
spontaneous symmetry breaking in this theory initiated by a quantum quench in has 
been rather extensively studied in the particle physics and cosmology literature (see for 
example Ref. |23j and references therein), where the principal motivation is to develop 
intuition for phase transitions in the early universe. 



■^Given the relatively short lifetimes of ultracold systems, it is also likely to be the only part of the 
story accessible to experiment for the foreseeable future. 

^We hope the meaning of classical is clear; superfluid hydrodynamics is a classical theory though 
superfluidity has a quantum origin. 
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There is a hierarchy of approximations of increasing complexity that have been apphed 
to this and related problems, reflecting the difficulty of extending the description of the 
behavior of the system to longer and longer times. As discussed above we will focus 
dynamics immediately following a quench, though we close this section by discussing the 
theoretical approaches that have been used to understand the late-time behavior. 

2. 3. Instability at the Gaussian level: emergence of classicality for one mode 

The very simplest approximation involves discarding completely the quartic term 
in Eq. Q, leaving a Gaussian model that may be solved exactly. Understanding the 
dynamics after a quantum quench in this model is still a useful exercise, however. Since 
the Hamiltonian separates into a sum of harmonic oscillator Hamiltonians for modes of 
each wavevector, we start by considering a single oscillator. 

H=Ip' + lu'x' (2) 

Let us suppose that at time t = the frequency is abruptly changes from a real to an 
imaginary value > — > < 0. The subsequent evolution of the operators x{t) and 
p{t) in the Heisenberg picture coincides with the classical solution 

xit) = :z:(0)coth(|c./|t) + g| sinh ^ ^ (^x{0) + (3) 

p{t) = p{0) coth (l^^li) + a;(0)|^/| sinh {\uf\t) > ^ (p(0) + \ujf\x{0)) . (4) 

Condensed matter physicists would be inclined to call the resulting time evolution the 
unfolding of a Bogoliubov transformation, while quantum opticians know this operation 
as squeezing. 

This solution may be used to evaluate any desired correlation function, using the 
initial correlation functions at time t = corresponding to the ground state with — uf. 
For example 

{x{tf) = coth^ i\u;f\t) {xiOf) + sinh^ {\u;ft\) (5) 



cc 
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c2|"/|t 



as i — oo. (7) 



The exponential growth is exactly what would be expected from the instability of the 
classical system, though of course {x{t)) = by symmetry. Quantum mechanics brings 
the novel feature that the instability is effective even in the absence of noise or some 
other perturbation to push the system away from the unstable maximum: the zero point 
motion in the initial state is enough. 

The quantum origin of the initial 'seed' notwithstanding, the description of the system 
becomes increasingly classical as time progresses. To understand why this is so, notice 
that the i — oo limits of x{t) and p{t) given in Eq. (§ are commuting. Now of course, 
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Figure 1: The Wigner distribution in Eq. JM with iiiit = 0.7, |a;j| = cji. The dashed Unes indicate the 
contours of the Hamihonian H = ^p^ + ^uj'j,x^. The initial distribution is axially symmetric but is 
'squeezed' along the classical trajectory with the passage of time. 



the evolution is unitary, so [x{t),p{t)] = i, but this lack of commutativity is increasingly 
irrelevant as time progresses, so that in the long time limit one can take p{t) — \u!f\'^x{t), 
with x(t) having a Gaussian distribution with the above variance. 

It's also useful to consider the Schrodinger picture counterpart of the evolution in 
Eq. ([3]). The initial Gaussian is annihilated by x{0) + ip{0)/(jJi, so the state at time t is 
likewise annihilated by x{t) + ip{t)/uji, corresponding to the wavefunction 



ipix,t) = A{t)exp ( ^^^^ tan{9 + i\ujf\t) 



-> A{t) cxp - 



lX^\ujf\ 



1 - 2e"'e 



-2\u,f\t 



(8) 



where tan^ = uji/\ujf\, and A{t) is a normalizing factor. The Wigner distribution of the 
associated density matrix p(xi,X2,t) = ip* {xi,t)'ip{x2,t) is 



f{p,x,t)= / dxi2r{x + Xi2/2,m{x~Xu/2,t)e-'P^- 



(X exp ~a{t)x' 



2 [p-PitM 



(9) 



a{t) 



as t — > 00 



with a{t) = 2|w/|sin26ie-2|"/l*, = |a;/| (l - 2cos26'e~2|"/l*). The form of this 

distribution is illustrated in Fig. [T] Because of the decay of a{t), the width of the x 
distribution is increasing exponentially, while at long times the p distribution is effectively 
S(p— \ujf\x). In other words, it is nonzero only on the phase space trajectory p = Iwjpa: 
of a classical particle that starts at p = a; = 0. 
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The exponential growth in the width of the x distribution also reflects the classical 
trajectory of the particle, so that the probabilistic character may be attributed to the 
initial state by writing 

x{t) = xoel"^!* 



where xq is a classical Gaussian random variable with variance (a;^^ — ^ 



(see 



0/ - 8 

Eq. ([5])). We may go further and attribute this randomness to fluctuations in the time 
taken to depart from the unstable maximum by defining a delay time t through |n] 



x = ±^J{xl)e^^f^^'-^l 
It is then a straightforward matter to find the distribution of r 

P(r) = y||c.;|e-l"/l-exp (^-^e-^l-/!-^ (10) 

known as the Gumhel distribution in probability theory, where it is often encountered in 
the statistics of extreme values or - as here - in survival analysis. 

After the emergence of the classical probability distribution, it is natural that the 
subsequent evolution be described by deterministic equations. The approach of combin- 
ing classical evolution with random initial conditions is often known as the 'truncated 
Wigner approximation' ^M]. The effect of finite temperature on the intial state can be 
included without difficulty. We will have more to say about the classical field approach 
in Section [2J1 



2-4-. Many modes: quench in the free field theory 

It's now straightforward to generalize the above results to the normal modes of N 
free quantum fields 



k 



-'(r) - ^ E XR^""' («L - ^ ^^l.■■■N (11) 



/Vol V V 2 

k 

The above expressions are for real fields: <?!>|(r) = 0j(r), 7r|(r) = 7rj(r). We can form 
the complex field that appears in Eq. ([ij from N — 2 real fields by taking $(r) = 
^ {(l)iir) + 102 (r)] and n(r) = ^ [7ri(r) - i7r2(r)]. We will stick with TV = 1 for the 
remainder of this subsection for clarity: the generalization to many fields is easily written 
down. 

Taking the dispersion relation to be o;^ = c^k^ +to^, we imagine changing the 'mass' 
term from positive to negative at time t = 0: mf > — > m?< 0, leading to imaginary 
frequencies in the range |k| < |m/|/c (wk appearing in Eq. (nTj) refers to the pre-quench 
value). We can use the result of Eq. ^ to write the second moment ((/)k(i)</'-k(i))j 
before Fourier transforming. Slightly neater expressions are obtained by considering the 
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Figure 2: Correlations between two points can only be affected after a quench once they lie in the forward 
light cone with apex halfway between them at the moment of the quench. 



'deep quench' limit in which rrii 3> I'ti/I, in which case we make take ajk,i = mi for 
|k| < \mf\/c, leading to the integral 



e^'' ('— ^Osinh^ (i^|m/|2-c2k2 
'2Try |m/|2 - c2k2 

The growth of the resulting correlations is determined by causality. To see this, note 



(0(M)0(r',t)) ^ J ^ ^-A-^j^^ ^. (12) 



that for |r — r'l > 2ct we may close the contour of integration in Eq. (12 1 in the upper 
half plane of |k|. There are no poles, so the integral vanishes. Correlations are absent in 
the initial state by virtue of the large rrii limit, and remain zero at time t while the the 
two points r and r' lie outside the forward light cone of a point halfway between them 
at i = (see Fig. [2]). 

Causality places similar (though less severe) restrictions on the propagation of corre- 
lations in non-rclativistic quantum mechanics, first studied by Lieb and Robinson [5]. 

2.5. Defect density in the Gaussian approximation 

The Gaussian approximation also allows us to address the formation of topological 
defects. Let us consider the case of vortices in a complex (or two component) field $(r, t), 
relevant to the two examples mentioned in Section |2.2| Usually one thinks of a vortex 
as an object with a well-defined core region, within which the amplitude of the order 
parameter goes to zero. Topologically, however, a vortex is a generic zero of a complex 
field, see Fig. [3j For Gaussian fields the density of such zeroes must be encoded in the 
second order correlation function, and is given by the Halperin-Liu-Mazenko formula 
[23 [26] 

n{t)^-^f"{r,t), (13) 

where f{r,t) = ($^(0)<i>(r))/(<i>^(0)$(0)) is the normalized correlation function. Apply- 
ing this result to Eq. ( [T2| ) (or rather its = 2 counterpart) gives immediately 

, ^ \fnf\ , ^ 

n(t) = ^^, asi^cx). (14) 




Figure 3: A complex Gaussian random field with correlator given by Eq. | |12| l at i = 2 (distance is 
measured in terms of the 'Compton wavelength' c/\mf\). Solid and dashed lines are the nodes of the 
real and imaginary parts of the field respectively, with vortices occurring at their intersection. Color 
indicates the phase of the field. 



It is interesting that even this simple model displays a decay in the density of defects 
as the order parameter coarsens, leading to vortex-antivortex annihilatior{^ Note that 
interpreting the quantum two-point function as a classical correlation function hinges 
upon our earlier discussion of the emergence of the classical picture. 

The spontaneous appearance of vortices in the transverse magnetization of a spinor 
condensate due to this mechanism were observed in the experiment of Ref. [15] that 
studied quenches in the spin-1 Bose condensate. The expected density of topological 
defects in this experiment was further discussed in Refs. [HI [171 HB] • As mentioned in 



Section 2.2 the other appearance of the relativistic 'At/)^' in ultracold physics is in the 
description of the superfluid to Mott insulator transition in the Bose-Hubbard model. 
Symmetry-breaking quenches from the Mott to superfluid phase were studied in a re- 
cent experiment |28j . though with no conclusive observation of superfluid vortices. A 
somewhat related observation of vortices formed during the Bose-Einstein condensation 
transition was reported in Ref. [29] . 

Taking a step back, we see that the spontaneous formation of topological defects upon 
passage through a phase transition is intimately connected to the finite propagation time 
of causal influences. Ordering is destined to proceed independently in regions sufficiently 
remote from each other, with topological defects a natural consequence. This general 
picture goes by the name of the Kihhle-Zurek mechanism (see e.g. Refs. ^30ii27j). Ref. [31] 
was an early test of this idea in a liquid crystal. 



'^This is the same behaviour as predicted by the Vinen equation for the decay of the density of vortex 
lines (see e.g. Ref. 127) ') but presumably for quite different reasons. 
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2.6. Kibble-Zurek scaling 

So far we have considered only instantaneous quenches. It is natural to ask what 
changes when the control parameter is allowed to vary at a finite rate. In this case a simple 
argument, due to Zurek, describes the overall scale of regions that order independently. 
Suppose that the control parameter S changes linearly with time, S = t/rg, with the 
phase transition occurring when (5 = 0. The correlation length at the transition diverges 
as ^ oc (5"'', while the characteristic relaxation time of the order parameter r is related 
to ^ by the dynamic critical exponent z, and thus displays critical slowing down r cx 

cx S^'"^. At some time before the transition, the relaxation time exceeds t^, so 
that further growth of the correlation length is impossible. This time scales with rg, 

the quench duration, as t* oc Tq'^'" , corresponding to a length scale oc Tq^'"' . This 
characteristic dependence on quench duration is known as Kibble-Zurek scaling. An 
extension of the Gaussian calculations of the previous section to the case of finite quench 
rate is consistent with this idea [3?. 

Considerable theoretical and numerical work has built upon the Kibble-Zurek idea [331 
1341 135[ 1361 137j , but some basic questions remain. Some derivations rest on the notion 
of "excitation" [35] or defect, whose meaning is clear in an integrable model such as 
the quantum Ising case but not in a generic non-integrable model; for at least some 
quantities such as entanglement, the behavior after a sweep in a non-integrable model 
is quite different (it increases more rapidly than in the integrable case). We are not 
aware of any experiment showing Kibble-Zurek scaling in an ultracold system. The 
most convincing tests to date were carried out in superconducting Josephson junctions 

m- 

2. 7. Putting on the brakes 

The Gaussian approximation used thus far is enlightening but naturally has its limits. 
Most obviously, the mode amplitudes in the band of unstable wavevectors satisfying 
c|k| < jmyl grow exponentially without limit. In addition, there is no scattering between 
different momentum states, the usual source of thermalization. In this section we briefly 
discuss approaches that allow these limitations to be overcome. 

It is clear that the inclusion of the quartic term in Eq. ([T]) will cause the instability to 
saturate. The Hartree-Fock approximation leads us to replace this term with a quadratic 
one 

^ -2(|$|2)2+4(|$|2)|$|2. (15) 

Note that there are no 'anomalous' terms ($2) — Q, reflecting the absence of symmetry 
breaking in the initial state. The resulting quadratic term can therefore be incorporated 
into the mass term to given an effective mass 

m2g = m2+4A(|<i>|2),„^„. (16) 

The right hand side is rUcs dependent, so this is a self-consistent equation. As in the 
stationary problem there are ultraviolet divergences to deal with. This method has been 
used a great deal to study field theories out of equilibrium: some early references are 
Refs. [ini 111] , while Ref. [42] is a recent study of the quench problem. 

The solution of the self-consistent Hartree dynamics saturates the instability, and 
may be formally justified in the limit of a large number N of field components (for 
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the complex field $(r) N = 2) with XN held fixed. But this is cold comfort, for this 
approximation ignores scattering between the different momentum modes, a sine qua 
non for thermalization. 

There are two contrasting pictures of the subsequent evolution. One could take the 
viewpoint of kinetic theory, using some approximation to close the hierarchy of equations 
of motion. The derivation of the Boltzmann equation is the most familiar example of 
this strategy. An enormous literature exists on the corresponding approaches to field 
theories. Ref. [15] provides a pedagogical introduction: it must be admitted that the 
high energy literature is considerably more sophisticated than the approaches typically 
encountered in treatments of these problems by condensed matter physicists. The pos- 
sibility of developing a kinetic description without the 'quasiparticle ansatzj^ could be 
important for the description of strongly correlated systems out of equilibrium. The 
second viewpoint focuses on the dynamics of the classical field distribution that emerges 
after the quench. Usually called the 'classical field method', this approach also has a long 
history, but has enjoyed particular popularity with the advent of Bose condensation in 
atomic gases, where systems are often dilute and well described by the Gross-Pitaevskii 
equation. 

Of course these approaches are related. Sticking with the example of the non- 
relativistic Bose gas (instead of the relativistic case Eq. ([!]), where the existence of particle 
and hole excitations complicates matters), one may derive a kinetic description starting 
from the classical equations of motion, in which the occupation number corresponds to 
the mean intensity of the field amplitude^ The very same problem has been studied 
extensively in the theory of nonlinear waves, where it is known as wave turbulence |44j . 
Closure of the equations is possible for weakly nonlinear, dispersive waves, in which case 
the amplitudes evolve to a joint Gaussian distribution |45l 146] . One expects, however, 
that the conditions of phase ordering (i.e. accumulation of particles in the low momen- 
tum states, in the case of Bose-Einstein condensation) are precisely those under which 
this approximation breaks down [TTl H5] . The introductory section of Ref. [JS] provides 
a useful summary of the difficulties involved in an analytical description of the kinetics 
of Bose condensation. In this case one would prefer a hybrid scheme in which the long 
wavelength components of the field ('condensate') are described by a classical equation 
of motion, with the short wavelength components ('normal fluid') being described by 
kinetic theory. In such a description the condensate must be subject to noise arising 
from stochastic interactions with the normal fluid. The derivation of such a 'two-fluid' 
picture also has a long history, see e.g. Refs. [SUl [SH [511 [531 \EM recent treatments. 

Finally, there is the matter of connecting such treatments to the extensive literature 
studying effective stochastic-dissipative models for equilibrium dynamic critical phenom- 
ena (notably the alphabetic classification of Hohenberg and Halperin [SSQ) and phase 
ordering [T5] . 



Essentially this amounts to a separation of scales in which the occupation numbers and spectral 
function of particles change slowly compared to the typical energy, see e.g. Section 4.2.2 of Ref. l:i 

*The resulting collision integrals only allow for stimulated, and not spontaneous scattering. The 
latter is negligible for highly occupied modes. 

^Ref. 1561 is a recent example showing that the classical field dynamics of the real scalar theory 
is consistent with the predictions of Model C, describing a system with non-conserved order parameter 
coupled to a conserved (energy) density, the latter being a feature of Hamiltonian dynamics. 
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The two incarnations of the 'At/)'*' Hamiltonian discussed in Section 2.2 are of recent 
vintage, but it is puzzhng that there has been relatively little fruitful interaction between 
theory and experiment p°]concerning the kinetics of Bose-Einstein condensation in atomic 
gases, a phenomenon occurring at the start of almost every experiment in ultracold 
physics. 

3. Integrability and some consequences 

As we noted in the introduction, ultracold physics provides a number of situations that 
come very close to realizing certain integrable models. It is a long-standing conjecture 
that finite dissipative transport coefficients - such as the mobility - can only emerge for 
systems that are non-integrable [51] |B2J . Atomic systems therefore form a natural 
testbed for this idea, in which one could hope to study transport as integrability is 
broken to varying degrees. Before continuing we should admit that even giving a precise 
formulation of the idea of quantum integrability is a notoriously thorny matter |64j . To 
be concrete, we adopt Sutherland's physical point of view that those integrable systems 
studied in condensed matter (usually by the Bethe ansatz, in one of its avatars) display 
non-diffractive scattering |65| . This means that when three or more particles undergo 
a collision, the only possible outcomes are those in which the set of outgoing momenta 
coincides with the set of incoming momenta, rather than taking on all possible values 
allowed by energy and momentum conservation. This statement is evidently equivalent to 
the existence of a much larger set of conserved quantities of the form ^ . /(fc,), for some 
arbitrary function of momentum. It is also equivalent to the factorization of many-body 
scattering amplitudes into successive two-body scattering event sj"^ 

As explained in Section [l] when this observation is translated into kinetic theory, the 
resulting lack of relaxation of the distribution function in a ID integrable system makes 
the absence of dissipative effects plausible. In this section we discuss this issue with the 
help of a simple concrete example. 

Let us consider an impurity moving in a spinless Fermi gas. When the impurity and 
gas masses are equal M — m this problem is integrable for any strength of the interaction 
^int — ~ ^) between the fermions at positions {xi} and the impurity at -Xp^ 

Low energy scattering events near the Fermi surface have momentum transfers close 
to a multiple of 2pp^ where pp is the Fermi momentum. Thus for an impurity of mass M 
there is a characteristic energy scale -Erccoii = '^Pp/M associated with a2pp momentum 
transfer. At temperatures ksT <C -Erecoii such processes are frozen out, and the motion 
of the impurity is determined by forward scattering. 

Consideration of the kinematics of these forward scattering events shows that pro- 
cesses involving a single fcrmion are also suppressed. Let us write the dispersion relation 
of the fermions as ^{p) = p^ /2m — ji for chemical potential /i, and that of the impurity 
as t{p) — p^ /2M . Scattering of a fermion with momentum p ^ ipp and low momentum 
transfer q leads to a change in energy of £,ip + q) ~ ^{p) ^ zLvpq. Since the corresponding 
change in energy of the impurity is on the order of g^/2M, energy conservation would 



If Notable exceptions are Refs. ISTl [58l l59l l60l l29l . 

^^For (5-function interactions this is more or less clear. 

^^This is a special case of the Gaudin-Yang gas, integrable for any number of 'gas' and 'impurity' 
particles [6611671 . 
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Figure 4: Diagrams contributing to the two-particle scattering amplitude in Eq. | |17| l. Solid lines denote 
the impurity, dashed lines the fermions of the gas. 



require q ~ 2Mvf = '^{M/m)pF, or energies ^ {M /mYE^ccoW- This is just a version of 
the familiar Landau criterion. 

These arguments suggest that the low temperature transport of the impurity is due to 
higher order processes. The simplest such process involves the scattering of two particles 
with momenta lying close to each fermi point. The lowest order amplitude for fermions 
with momenta ki and k2 to scatter to ki +qi and ^2 +92, while the impurity momentum 
goes from K io K - qi~ q2 (see Fig. H is T^'^^S{E, - Ef) with [32] 



7-(2) 



1 



1 



1 1 



+ 92 



Cfci 



+91 



e/i'-fei-gi+fc2 



ex 



(17) 



At low gi, q2 the first and third terms give rise to singular behavior arising forward 



scattering processes at second order T^^^ ~^ * (r) vlqiq2 
momentum relaxation rate of the impurity, given by 



. Despite this singularity the 



27r 
hMT 



J2 («i + 92)'|r(2)p(5(i?. - Ef)nk,n,^ (1 - n^^+^J (1 - n^.+^J 



ki,k2,qi,q2 

{jik is the Fermi distribution) is finite and vanishes as T^, leading to a low temperature 
mobility = Tnwm/M cx T^'^ (SHI [ZO] • This calculation, together with the observation 
that even an almost opaque impurity will appear transparent at low temperatures as the 
backscattering processes are suppressed, offers a qualitative picture of behavior of the 
mobility from high to low temperatures. |70j 

However, the above matrix element in fact vanishes in the limit Af = to of equal 
masses, due to interference between the terms. This corresponds to the absence of 
diffractive scattering mentioned at the beginning of this section, and coincides with the 
integrable limit. A related observation was made by Gangardt and Kamenev |71j in the 
case of the Bose gas, where the process analogous to the creation of two soft particle-hole 
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pairs is the creation of two phononj^ 

But what about the two-particle scattering processes? Though they may be expo- 
nentially suppressed at low energies, they can certainly change the momentum of the 
impurity. It is not immediately clear from this picture why ballistic motion of the impu- 
rity is guaranteed when M = m. To understand why this is so, we use a variant of the 
Kohn formula for the Drude weight [72l [62] to express the singular part of the mobility 
Re/i(a;) = fjDS{io) H as 

n 

where the sum is over the many-body eigenenergies En and Pn are the corresponding 
Boltzmann weights. $ is a flux felt only by the impurity particle. For a free particle 
D = 1. Now we invoke another feature of integrable systems: the occurrence of level 
crossings in the spectrum with the variation of a parameter that preserves integrability, 
of which $ is an example. Because of the absence of avoided crossings, it is plausible that 
the average curvature is nonzero in the thermodynamic limit. When AI ^ m, avoided 
crossings occur and D ^ may result. Ref. studied a lattice version of this model via 
exact diagonalization and the Bethe ansatz in the integrable case, finding some support 
for this conjecture. 

We are left to ponder the connection between nonzero D and non-diffractive scatter- 
ing. Some insight is gained by considering the curvature of the free energy with respect 
to the flux [73^ 



^ ^ hm i rdt((P(i)P(0))-(P)2) 



(19) 

ML^ {MLY T^oo T Jo 
d^F/d^^ = in the thermodynamic limit, so 

A iin,((p(i)p(o))_(P)2). (20) 

IVl t— >oo 

Thus a non-zero D corresponds to non-decaying momentum correlations. By Khinchin's 
theorem [73], this amounts to a failure of ergodicity for the variable P. The relation 
between this behavior and the conservation laws is provided by the Mazur inequality [75] 

lim mt)pm {p?) > E irnfs^ (21) 

where Qa are the conserved quantities, orthogonal to each other under the inner product 



provided by (•■■). Eq. (21 ) is in fact an equality if all conserved quantities are included 



^''To see the same thing happening in simple perturbation theory, one has to consider the integrable 
point of the bosons plus impurity problem where the boson-boson interaction strength is the same as 
the boson-impurity interaction, and both the bosons and the impurity have the same mass. 
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Doublon states at V=l , L=10 




Figure 5; Spectral signature of doublon decay in the ID extended Hubbard model, with on-site interaction 
U /t, and nearest neighbor interaction V/t breaking the integrability of the model. The steep diagonal 
lines correspond to states of the system with a single doublon; the shallow lines to states with no doublon. 
The mixing between the two types of states that results in avoided crossings is the spectral signature of 
doublon decay that would be absent at y = 0. 



[76] , The resulting connection between ballistic transport, the (thermally averaged) level 
curvature, and the presence of conserved quantities has been explored in a number of 
works: the introductory sections of Ref. [77j provide a detailed summary. 

We close this section by pointing out that the above ideas are of more general appli- 
cability. Consider the example of doublon decay, the decay of doubly occupied states in 
the fermion Hubbard model, already observed in ultracold gases [78]. In the ID Hubbard 
model 



Hubbard 



= -t 



u 



(22) 



this process is forbidden by integrability, but it can occur when integrability is broken. 
The doublon number operator is Nd = X^i '^t,*"'^,*' where Us^i denotes the number of 
spin s particles on site i. But since this is just the on-site interaction term the above 
reasoning implies 



^.Pn^ < 13 lim {{Nd{t)Ndm 



dm 



(23) 



where we now have an inequality because the 'susceptibility' —d^F/dU^ should be posi- 
tive for stability. Doublon decay is thus associated with a vanishing of the averaged level 
curvature with respect to U, see Fig. [5] 
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4. Disordered systems and many-body localization 

A fertile area of condensed matter research has been the physics of systems with 
"quenched" disorder: there are random quantities in the Hamiltonian that are not in- 
cluded in the thermal average. One may be interested in either the mean values of 
physical quantities after some non-thermal average over the disorder realizations, or in 
typical values that are the most likely to be observed even if rare realizations dominate 
the average. The possibility of retaining quantum coherence for long times in ultracold 
atomic systems may allow experimental answers to some basic questions in the theory 
of localization in quantum disordered systems. The 'Anderson localization' of a single 
particle moving quantum-mechanically in a random potential is in many cases well under- 
stood. Remarkably a particle moving in one or two spatial dimensions has only localized 
energy eigenstates (i.e., the wavefunction falls off exponentially at spatial infinity) even 
for a weakly disordered potential; this is in contrast to the case of three dimensions, 
where a "mobility edge" separates localized wavefunctions at low energies from extended 
wavefunctions at high energies. 

The effect of interparticle interactions on this description has been a subject of active 
research for many years but is still not completely understood |79] . For experiments 
on electrons in solids, the standard approach for describing experiments is to assume 
that an electron loses its phase coherence after some temperature-dependent dephasing 
time. This dephasing time provides a cutoff on the localization physics, as localization is 
a consequence of phase coherence. Perturbative treatments of interparticle interactions 
have had success in describing some experiments, but at least one basic question remains 
unanswered. 

The dephasing time is an example of interaction with a "bath" , similar to those 
discussed earlier: the electron system in a disordered metal is treated as open rather 
than Hamiltonian. It is natural to ask what happens in an interacting particle system at 
nonzero temperature when there is no external bath. (Since there is no bath, the meaning 
of "temperature" is simply a weight on initial conditions for Hamiltonian dynamics, as 
in the Kubo formula.) From our discussion of equilibration, one might think that one 
part of the system can see the rest of the system as such a bath, at least approximately, 
so that interactions alone lead to an apparent dephasing time. 

Recent theoretical work suggests that there can be a localization transition in a 
closed (i.e., Hamiltonian) interacting system at some finite temperature Tf,, when all 
single-particle states are localized [TT]. Above Tc the system acts as its own bath and 
appears dephased on long length or time scales, while below Tc the localization persists. 
Actually the problem can be simplified somewhat by switching from continuum 'particle' 
systems with an unbounded entropy per unit volume at infinite temperature to lattice 
systems with finite maximum entropy per site [80] . Then the many-body localization 
transition exists as a function of parameters even at infinite temperature, and we will 
discuss an example believed to exhibit such a transition. 

We start from an XX spin chain with random Zeeman magnetic fields: 

Ho = -JY^ (Sf + SfSf^^ - J2 h^St. (24) 

i i 

Here the fields hi are drawn independently from the interval [—A, A]. The eigenstates 
of this model are equivalent through the Jordan- Wigner transformation to Slater deter- 
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minants of free fermions with nearest-neighbor hopping and random on-site potentials. 
Since every single-fermion eigenstate is locaHzed by those potentials as long as A > 0, 
we expect that the dynamics of this spin Hamiltonian are localized as well. The response 
of the system to a local perturbation is confined for all times with exponential accuracy 
to a radius set by the single-particle localization length. 

Now if instead we had started with Heisenberg rather than XX couplings between 
the spins, i.e., added a term —J^^S^Sfj^-^, then the model, at infinite temperature, 
is thought to have a dynamical transition as a function of the dimensionless disorder 
strength A/ J [5DJ[5T]. This transition should be manifest in several physical quantities, 
including level statistics and dynamical correlation functions. Level statistics should 
transition from Wigner-Dyson statistics and level repulsion for small A/ J to the Poisso- 
nian distribution characteristic of a disordered system (since spatial regions separated by 
the localization length are effectively independent) for large A/ J. Correlation functions 
of the conserved quantity Sz (equivalent to particle number) should show, as originally 
emphasized by Anderson (10) . a transition in the long-time behavior between particles 
escaping to infinity or remaining in a bounded region. 

A system in the localized phase may fail to come into mutual thermal equilibrium 
among its constituent pieces. To understand why this is so, consider dividing the above 
spin chain into two equal halves, A and B. Though the disorder is not the same on each 
side, we expect the two halves to be macroscopically identical. Conventional thermody- 
namic wisdom would have that in the long time limit, the expectation value of total spin 
S^jg = Yli^A/B t^*^ same on both sides. 

Let us take the initial state of the system to be 

1*) = \Ma)a ® \Mtot - Ma)b (25) 

where are states of the A and B subsystems with S\^g — Mj^/g. Evidently 

there are many such states. In terms of the eigenstates of the whole system, the long 
time average of AiS^ — S% — Sg is then 

(A5^ = 5](n|A5^|n)|(n|*)p. (26) 

n 

In this sum the first factor refiects the distribution of A5^ in the eigenstate while 
the second factor reflects the overlap of this state with the initial state. 



Let us consider the uncoupled system, where the appropriate hopping term in Eq. ( 24 ) 
that connects the two regions is set to zero. The eigenstates \n) are then eigenstates of 
A5^, and no relaxation occurs. The vast majority of states have AiS^ eigenvalue AM ^ 0, 
and the distribution of AM generally has the form Ps(AM) = exp [-iVs( AM/TV)]: 
deviations from the equilibrium value are exponentially rare. The s subscript is to remind 
us that the function s{x) is just the entropy (per site). 

When a small coupling J is switched on, states with different AAf are mixed to- 
gether, but we expect them to retain their character, meaning that the histogram of AAf 
probabilities associated with a particular state will be strongly peaked around the value 
that it had with certainty for J = 0. Similarly, the probabilities | p will be largest 

for those states peaked around AMq = 2M^ — Mtot, the initial value. Let us group the 
eigenstates in sets Mam labelled by the value of AAf that they had when J = 0. Then 
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defining, P^{^M) = E„ga/-am I ^q- ^ has the form 



(A5^) - ^ AMP,{AM)P^{AM). (27) 
aj\/ 

Since states in Nam only overlap with |^') at the \AM — AMo\^^ order of perturbation 
theory in J, it is reasonable that P^;{AM) is exponentially small is this deviation. Thus 



Eq. (27) is seen to involve a competition between two exponential factors that in the 



thermodynamic limit can pick out a value of (AS^) different from zero, contradicting the 
assumption that the two parts of the system come to mutual equilibriurrj^ 

The crucial part of the picture is that the states 'retain their character' at J ^ 0, 
and this is just what is implied by the recent work discussed above. When there is no 
longitudinal part to the interaction, this is just Anderson localization of noninteracting 
Jordan-Wigner fermions. What is new is the claim that the same remains true in the 
interacting problem over some range of J. 

However, the location of the transition and its critical properties have been difficult 
to obtain even numerically. A search based on level statistics was inconclusive [80] . 
while studying correlation functions suggests a phase transition [STl In both cases ex- 
act diagonalization of the Hamiltonian was used, which strictly limits accessible system 
sizes. The phase transition in the correlation functions was interpreted in terms of an 
"infinite-disorder" critical point, analogous to the critical points in ground-state quan- 
tum phase transitions in disordered systems accessible via the real-space renormalization 
group (RSRG) [83l|84]. A approach to many-body localization based on numerical solu- 
tion of RSRG equations predicts a sharp transition and a specific location [85j , but it is 
currently difficult to compare these predictions to numerics on the microscopic model. 

Specific proposals for realizing a many-body localization transition with ultracold 
atoms are discussed in Ref. [5^. Another way of searching for the putative many-body 
localization transition, at least numerically, involves the entanglement between different 
spatial regions. In the disordered phase, one would naively expect the entanglement 
entropy, defined in a moment, to saturate, as indeed happens in the non-interacting case, 
while in the extended phase entanglement would spread throughout the system. We 
discuss an interesting behavior in the dynamics of entanglement in the disordered phase 
below, after explaining some general features of entanglement in many-body systems. 



5. Entanglement dynamics 

One of the most significant potential advantages of interacting many-atom systems is 
the ability to maintain quantum coherence for relatively long times than in the electron 
subsystem of solids. Considerable effort has been devoted in recent years to under- 
standing how entanglement, one of the basic notions of quantum information, behaves 
in many-particle systems, both for intrinsic interest and for potential application to 
quantum computing. Since atomic systems at the moment may offer the best hope for 
observing many-particle entanglement in a quantitative way, and since entanglement un- 
derlies many of the questions touched on in preceding sections, we now review the basic 
notions of entanglement and its dynamics in model atomic Hamiltonians. 



The importance of rare states in the absence of thermalization was emphasized in Ref. |82| 

18 



5.1. Basics of entanglement in ground states of many-body systems 

In our discussion of entanglement, for simplicity we will limit ourselves to the case 
of a pure state of a bipartite system AB. While the Hilbert space of the full system is 
spanned by product states (^\4>j) of basis states of A and B, there are superpositions 
of the basis states that cannot be factored into any pure states of A and B. A familiar 
example is the singlet state of two spins: \'^ab) — ^ (I t^) <8i | Ib) — \ ts) ® | Ia)) ■ An 
entangled state is simply one that is not a product state. The entanglement entropy, 
which has been the quantity most studied in the context of many-particle systems, is 
defined as the von Neumann entropy of the reduced density matrix for either subsystem, 

S = -TrpA log Pa = -Trps log , (28) 

where pA (pb) is the reduced density matrix of subsystem A{B) obtained by tracing over 
degrees of freedom in the other subsystem. 

The basic behavior of entanglement entropy in ground states of standard atomic or 
condensed matter is generally well understood, especially in one dimension. We will not 
discuss ground-state entanglement entropy in detail here but simply state a few results 
before referring the reader to recent in-depth reviews. In gapped systems, the general 
expectation, which can be proved in some cases, is that entanglement entropy satisfies 
an "area law", i.e., scales as the volume of the AB boundary. For regular geometries, 
the scaling is thus as L'^~^ for a subregion A of linear size L cut out of an infinite d- 
dimensional system. Topological phases can have subleading terms in the entanglement 
entropy that are of order unity (i.e., scaling as L'^), in addition to the area law, and probe 
aspects of the topological order [87l [88] . 

For critical systems, the behavior of ground-state entanglement entropy is more com- 
plicated. In one dimension, for quantum critical points described by 2D conformal field 
theories, the entanglement entropy diverges logarithmically, with a universal coefficient 
determined by the central charge [HHl [HHl IE] : 

Q 

SKi-\ogL/a as L -> cx) (29) 

o 

for a block of length L cut from an infinite chain, where c is the central charge and a a 
short distance cutoff. There are similar logarithmic divergences at infinite-randomness 
critical points in one dimension [84], although differences between the pure and random 
cases |92j become manifest in the entanglement spectrum |93| . In higher dimensions, 
critical points with a Fermi surface can violate the area law [MJ [US] , while others can 
obey the area law but also have subleading corrections of interest [S5] . 

To date it has been difficult to observe entanglement directly except in small sys- 
tems |97j . Since many-particle systems with local Hamiltonians are expected to have en- 
tanglement described by the area law, or possibly larger by a logarithmic factor, thermal 
entropy (which scales as volume L'^ under the circumstances above) will rapidly domi- 
nate except at the lowest temperatures. However, there are recent proposals (HSJ IHS] for 
how relatively standard measurement processes could probe many-particle entanglement. 
Entanglement also determines the ability of certain numerical methods based on matrix 
product states and their generalizations to capture accurately the ground state of a given 
Hamiltonian. 
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5. 2. Entanglement dynamics of low-lying states 

The growth of entanglement after a quench or sweep of parameters in the Hamihonian 
has been studied for one-dimensional systems using similar methods to those studied 
above in Sections [2j For a quench in a conformally invariant system, one can distinguish 
two simple scenarios: a "global quench" in which the parameters are changed instantly 
everywhere in space, versus a "local quench" in which the parameters are only changed 
at one point [H]. In both cases, entanglement propagates along "light cones"; this can 
be understood intuitively by noting that conformal invariance implies that all excitations 
propagate with a single velocity. For a global quench, this leads to entanglement growing 
linearly in time after the initial 

The same argument can be applied to understand the effects of sweeping across a 
conformally invariant critical point. We focus on the case of a translation-invariant 
problem, which is similar to the global quench as far as how entanglement propagates: 
after the sweep, if the Hamiltonian is held constant then entanglement entropy increases 
linearly in time. The rate of this linear increase, however, is now determined in simple 
models by the number of excitations that were created above the ground state when 



the energy gap was small (we make the same assumptions here as in Section 2.6). This 
number is determined as a function of the sweep rate F by the same scaling arguments 
as presented above |100j . Right after the sweep, the entanglement scales as a logarithm 
of the sweep rate with a universal coefficient combining central charge c and correlation 
length exponent v |101l 157] ; for the half-chain entropy, 

^ = - r + const, (30) 

In more complicated models with interactions between the "excitations" in the final 
gapped region of the phase diagram, the behavior of entanglement during a sweep is more 
complex. Without conformal invariance, there is no general rule known for the growth 
of entanglement with time, and indeed entanglement seems to grow quite rapidly once 
excitations begin to interact [ST- (Numerical studies in this regime become technically 
challenging as the entanglement quickly saturates the amount that can be captured with 
matrix product states of a given matrix size.) 

The growth of entanglement can be understood as a consequence of the "entanglement 
thermalization hypothesis" [1021 1103j : if the interactions lead to apparent equilibration 
for local measurements with an effective temperature determined by the (constant) en- 
ergy density, then the entanglement entropy of a subregion must reproduce the thermal 
entropy at that temperature, since this entropy can in principle be measured. But the 
thermal entropy will satisfy a volume law rather than an area law, and hence will be 
parametrically larger for a large subregion than the initial entanglement entropy, as the 
system began close to a ground state satisfying the area law (possibly with logarithmic 
corrections) . 

In closing, we note that the putative many-body localized phase discussed above also 
has a signature in the dynamics of entanglement that is not yet well understood. In a 
localized phase of a non-interacting system, an initial product state has entanglement 
that grows in time but saturates once the individual particles have diffused over a length 
scale set by their localization length. In an extended phase, one expects entanglement 
to increase more rapidly, presumably as some power law of time depending on whether 
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the transport is diffusive or ballistic. In interacting disordered systems, numerics sug- 
gest [1041 1105) a regime where the entanglement entropy increases logarithmically in time 
without limit. Whether this is generic behavior and indicates that the phase is not truly 
localized, or instead is a qualitative difference between localized phases with and without 
interactions, are questions for future theoretical work. 
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